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NOMENCLATURE 
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width of beam and piezo-actuator, m 

electric charge constant of piezo-actuator, m/v 

distance to neutral axis of composite beam measured from its 

lower edge, m 

Young* a modulus of elasticity of piezo-actuator and beam 
respectively, N/m 2 

Young* s modulus of elasticity of element i of the beam, N/m 2 

vector of modal forces and moments, N or Nm 

vector of external forces and moments, N or Nm 

area moment of inertia of actuator, bonding layer and beam about 

the neutral axis of the oomposlte beam respectively, m* 

area moment of inertia of the element i, m* 

mass moment of inertial of the composite beam at code i, kg-m 2 

stiffness matrix of the element i 

overall stiffness matrix of beam-actuator system 

length of element i, m 

mass of the composite beam at node 1, kg 

mass matrix of beam-actuator system 

external moment acting on node of beam, Nm 

piezo-electric moment generated by piezo-film, Nm 

number of elements of the beam 

thickness of piezo-actuator, bonding layer and beam respectively, 

m 

modal coordinates of the flexible system 

voltage applied across the piezo-electric film, volts 



external force acting on the 1 th element of the beam, N 
y^ the linear translation of node 1, m 

w b,f,s boDd » and beam maBS ^ er unlt length respectively, kg/m 

Greek Letters 
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•• 

*i 
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factors 1 or zero if actuator is bonded to element i or not 
respectively 

deflection of node 1, n or rad. 

deflection vector of all nodes of the beam, m or rad. 

acceleration vector of all nodes of the beam, m/s 2 or rad/s 2 

piezo-electric strain in piezo-actuator, m/m 

angular deflection of node i, rad. 

diagonal matrix of the eigenvalues of the system 

bending stresses in piezo-actuator, bonding layer and beam 

respectively, N/m 2 

piezo-electric stress in actuator, N/m 

modal shape matrix of the eigenvectors of the flexible system 



This study deals with the utilization of piezo-electric actuators in 
controlling the structural vibrations of flexible beans. 

A Modified Independent Modal Space Control (MIMSC) method is devised 
to select the optimal location, control gains and excitation voltage of 
the piezo-electrio actuators in a way that would minimize the amplitudes 
of vibrations of beams to which these aotuators are bonded, as well as 
the input oontrol energy necessary to suppress these vibrations. 

The presented method accounts for the effects that the 
piezo-electric actuators and the bonding layers have on changing the 
elastic and inertial properties of the flexible beams. 

Numerical examples are presented to illustrate the application of 
the MIMSC method and to demonstrate the effeot of the physical and 
geometrical properties of the bonding layer on the dynamic performance of 
the actively controlled beams. 

The obtained results emphasize the importance of the devised method 
in designing more realistic active control systems for flexible beams, in 
particular, and large flexible structures in general. 


UTRODOCTIOI 


Active vibration control systems are becoming essential and viable 
means for minimizing the vibrations of large flexible structures which 
are intended to provide stable bases for precision pointing in space. 

Distinct among the presently available active control systems are 
those that rely in their operation on piezo-electric actuators. Such 
systems have proven to be experimentally effeotive in controlling the 
vibrations of simple structural elements such as rectangular beams [1-2] 
and hollow cylindrical masts [33* The effectiveness of these systems is 
coupled also with the light weight, high force and low power consumption 
capabilities of the piezo-electric actuators [4-8]. These features 
rendered this olass of actuators to be an attractive candidate for 
controlling structural vibrations. 

The present state-of-the-art of this type of actuators has been 
limited to the analysis and testing of their characteristics [9-11] as 
influenced by their geometrical or operational conditions. A recent 
attempt has been made by Baz [12], to select their optimal geometrical 
parameters and location which are best suited for a particular structure 
subjected to known static loading conditions. The developed synthesis 
procedure has proven to be essential to the successful integration of the 
actuators into the structure in order to minimize its static deformation. 

However, no effort has been done to optimize the control of the 
vibrations of a multi-mode flexible system using a small number of 
optimally placed piezo-actuators through the development of efficient and 
realistic control algorithm that ensures minimal amplitudes of 
oscillation and input control energy. 


It is, therefore , the purpose of this study is to devise suoh an 
optimal oontrol method that is based on a modification of the well known 
Independent Modal Space Control(IMSC) method [13-15]. The devised method 
accounts for the spillover from the controlled modes into the 
uncontrolled modes due to the use of fewer actuators than the modeled 
modes. The developed method Incorporates also an optimal placement 
procedure that will enable the selection of the optimal location of the 
piezo-electric actuators in the struoture in order to guarantee a balance 
between the suppression of the amplitudes of vibration and the input 
control energy [16]. 

Furthermore, the developed procedure considers the effect that the 
piezo-eleotrle actuators and bonding layers have on ohanging the elastic 
and Inertial properties of the struoture to which they are bonded to. 
Suoh ohanges result in modifying the normal modes of the structure one 
way or another depending on the physical and geometrical properties of 
the actuators and the bonding layers. 

In this study, the emphasis will be placed on piezo-electric 
actuator-beam systems to illustrate the devised oontrol strategy. But, 
without any further modifications, the developed method can be readily 
applied to large structures. 

THE PIEZO-ELECTRIC ACTUATOR -BEAM STSTEH 
A. GENERAL LAYOUT 

Figure (1) shows a general layout of a flexible beam (A) whose 
deflection is to be controlled by piezo-electric actuator (B) bonded to 
the beam by bonding layer (C). The beam, under consideration, can 




generally be made of several steps which are not neoessarily of the sane 
thickness or the same material. The interfaoial nodes between the 
different steps can be subjected to external forces, moments or both. 
Further, the degrees of freedom of any node can be limited to linear 
translations, angular rotations or restrained completely depending on the 
nature of support at the node under consideration. 

In this study, the beam is assumed to have rectangular oross section 
of oonstant width b and that its transverse defleotion is due to the 
flexural action of the external forces and moments. 

In Figure (1), the piezo-electric actuator B is shown bonded to the 
element i of the flexible beam to form a composite beam. When an electric 
field is applied across the film, then it will expand if the field is, 
for example, along the polarization axis of the film and will contract if 
the two were out of phase. The expansion or contraction of the film 
relative to the beam, by virtue of the piezo-electric effect, creates 
longitudinal bending stresses in the composite beam which tend to bend 
the beam in a manner very similar to a bimetallic thermostat. 

With proper selection, placement and control of the actuator, it 
would be possible to generate enough piezo-electric bending stresses to 
counter balance the effect of the exciting forces and moments acting on 
the beam in a way that minimizes its structural vibrations. 

B. MODEL OF AM ACTUATOR-BEAM ELEMEHT 

Figure (2) shows a schematic drawing of a piezo-film A bonded to an 
element B of the flexible beam by bonding layer C. 

If a voltage v is applied across the film, a piezo-electric strain 
£f is introduced in the film and can be computed from : 




where d is the electric charge constant of the film, m/v 

t-j is the thickness of the piezo-electric actuator, m 
This strain results in a longitudinal stress Of given by : 

Of * E](d/tj)*v (2) 

where Ej is the Young’s modulus of elasticity of the film, N/m^ 

This, in turn generates a bending moment Mf, around the neutral axis of 
the composite beam, given by : 


Hf « 


-(t-i+to+ta-D) 

/ o f tb*y)dy 
-(t^t^-D) 


(3) 


where t 2 is the thickness of the bonding layer, m 

tj is the thickness of the beam, m 

b is the width of the beam, the bonding layer and the 

piezo-film, m 

In equation (3), D is the distance of the neutral axis from the lover 
edge of the beam which can be determined by considering the force balance 
in the longitudinal direction X of the beam, or : 
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-(to+to-D) - (t o -D ) D 

E«b f ydy + E 2 b J ydy + E 3 b / ydy s 0 
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( 4 ) 


( 5 ) 


where 


E 2 is Young’s modulus of elasticity of the bonding layer 
Ej is Young’s modulus of elasticity of the beam 


Equation (5) yields tbe following expression for D : 

E^t ^ 2 ^+E jt j^+2E ^ t ^ (t 2 +t 3 )+2E 2 t 2 t 3 

2(E^t^+E2t2+Ejtj) 


( 6 ) 


Equation (2), (3) and (6) can be combined to determine tbe bending moment 
Mf generated by the piezo-film on the composite beam as follows : 

d«b»E 1 «v»(E 2 t 1 t 2 +E 3 t 1 t 3 +E 2 t 2 2 +2E 3 t 2 t 3 +E3t 3 2 ) 
2»(E 1 t 1 +E 2 t 2 +E 3 t 3 ) 


Mf * 


(7) 


For this composite beam, it can be easily shown [12] that it has a 
flexural rigidity (E^I^) given by : 


Ejl^ t E^Ii + E0I9 + Bol 




2 A 2 * 6 3 A 3 


( 8 ) 


where Ij , I 2 and I 3 are tbe area moments of inertia of the film, the 
bonding layer and the beam about the neutral axis respectively. 

Let us now assume that the composite beam, shown in Figure (2), 
extends a length 1 A between two nodes (i) and (i+1). Further, it is 
assumed that the external forces and V 1+1 as well as the external 

moments M e ^ and M e ^ are acting on the beam at nodes 1 and i+1 
respectively. Then, the resulting linear and angular deformations of the 
beam y^ and 0 ^ as well as y^ and 0j^ at the codes i and i+1, 
respectively, can be related to the loads acting on the element as 
follows [173 s 
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(9) 



Equation (9) can be rewritten as : 


Fi = (10) 

where F^ la the resultant forces and moments vector acting on the beam 
element 1, N 

is the stiffness matrix of the composite beam element 1, N/m 
®i is the deflection vector of the nodes bounding the beam 
element , m 

Equation (9) constitutes the basic finite element model that relates 
the external loads (? and Mg) and piezo-electric moments (Mf) to the 
deflections (7 and 6 ) of the element as a funotion of its elastio and 
inertial parameters. 

The equation can be equally used for any element of the beam whether 
it has a piezo-film bonded to it or not. In the latter case, Mf is set to 
zero and flexural rigidity Ejl^ is set equal to that of the flexible 
beam element under consideration. 

The force-displacement characteristics of the individual elements of 
the beam-actuator system, as given for element i by equation (9)» are 
combined to determine the overall stiffness E of the beam system as shown 
in References [18-19) for example. 

The inertial properties of the composite actuator-beam system are 
determined using the lumped mass method where the mass and rotational 
inertial of each element is distributed among the nodes bounding the 
element [20]. 

Therefore, the diagonal mass matrix M (2n*2n) for the actuator-beam 
system, shown in Figure(l), can be written as : 
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where o 1 * [(w s +7 1 w f +7 1 w b ) f L 1 ]/2 

J 1 * t<V 7 1 w 

■l « t(w 8 +7 i . 1 v f +7 i . 1 w b )»L i _ 1 +(w b +7 i W f +7 i w b )«L i ]/2 
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“N+1 * t( w s +7 N w f +7 N w b )#L N^ /2 

J N+1 s [(w b +7 N w f +7 N w b )»L N 3 3 / 12 

v b,f,s s bond * film and beam mass per unit length respectively, 
kg/m. 

7^ * 1 if an actuator is bonded to beam element i. 

0 if not. 

The stiffness and mass matrioes K and M are used to define the 
dynamic equations of motion of the actuator-beam system. 

The equations of motion of the actuator-beam system can be written 


as follows : 


M£ + K5 = F 


( 12 ) 


nodes 


^(lj) is the nodal shape at node i and looatlon lj. 
The above equation can be rewritten as : 



If only C nodes are controlled with equal number of control forces Fq, 
then F R =0 and equation (6) reduces to : 


f c = B CC F C (18) 
and 

f R a b RC f C * 19) 

In the IHSC net hod, it is assumed that the control forces F c will 
not contribute to the excitation of the residual higher order modes. 
Accordingly, it was assumed that there is no control spillover from the 
controlled modes into the uncontrolled modes. Mathematically, this means 
that the IMSC method assumes that f R =0. This of course can only be true 
if the number of controlled modes is very large compared to the number of 
residual modes or when the residual modes are at much higher frequency 
band than the controlled modes. If these two conditions are not 
satisfied, then there will be considerable interaction between the 
controlled and residual modes. 

The MIMSC method considers such interaction by calculating the 
optimal modal control forces Cf q] using the IMSC close form solution of 
the Riccati Equation such that the control force f A of the i th mode, as 
given by [7 ) * is : 


fi 8 -(g-jcajUj ♦ g 2 Ui>/R 


( 20 ) 


where R is a factor that weighs the Importance of minimizing the 

vibration with respeot to the control forces, 
is the resonant frequency at the i tl1 r.-.mal mode. 
u if are the modal displacement and velocity respectively, 
g^, 62 are the modal position and velocity feedback gains given by 
[133 as : 

g 1 = -c^R +((w i R) 2 ♦w 1 2 R) 1/2 (21) 

g 2 - (2Ru^(-c»^R + ((wjR) 2 4^ 2 R)1/ 2 ) + Wi 2 R)1/ 2 (22) 

Accordingly, the displacement u A and velocity u A at the i* h mode can 
be feedback and used along with equations (20), (21) and (22) to 
determine the modal oontrol foroe f±. 

Once these forces are calculated, equation (7) is solved to give the. 
physloally applied control forces Fq as follows : 

F C * B CC**^ * (23) 

Then equation (8) is used to calculate the modal forces f R that 
would excite the residual modes which are generated by the spillover from 
the controlled modes. Definitely these f^ are not equal to zero as 
originally assumed in the IMSC method. 

Equations (14) can then be integrated with respect to the time to 
determine the modal displacements (u A ) and velocities (uj) which can, in 
turn, be used again to compute the modal forces f and so on. 

From the modal displacements and velocities, the physical state (6) 
of the flexible system can be determined from equation (13). A 
relationship can therefore be established between the physical state of 
the system and the physical control forces Fq applied to it. 


The MIMSC method Incorporates also an extremely Important feature 
vhioh is based on a "TIME SHARING" strategy of a small number of 
actuators in the modal space to control large number of modes. 

In this strategy, the modes of vibrations of the flexible system are 
ranked according to their modal energy level. If C actuators are to be 
used, then these actuators will be dedicated, at any Instant of time, to 
control the C modes that have the highest modal energy. In this way, the 
aotuators will first attenuate the modal energy of the oontrolled modes. 
During that time the control spillover will excite the uncontrolled 
modes. When the modal energy of the uncontrolled modes starts exceeding 
that of the controlled modes, the actuators are svitohed to oontrol these 
high energy modes in order to damp out their vibrations. Such time 
sharing of the actuators between the modes will eventually bring all 
these modes under oontrol. 

Figure (3) outlines a flowchart of the MIMSC method indicating the 
main steps of optimal placement and time sharing of the actuators as 
well as the consideration of the spillover between the controlled and 
residual modes. 


NUMERICAL EXAMPLES 


The application of the MIMSC method is illustrated by considering a 
three-element straight steel cantilever beam, shown in Figure (4), which 
is 0.0125m wide, 0.0021m thick and 0.15m long. A transverse impulsive 
loading of magnitude of 0.1N is applied at the free end of the beam for 
1.0ms. 

The effects of varying the thickness and material of the bonding 




Figure (3) - Flowchart of the 


MIMSC algoritha 
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Figure (t) - Schematic drawing of a three-element oantilever beam 



layer as veil as the location of the piezo-actuator are considered. 

The physioal properties of the considered bonding layers and a 
typical piezo-actuator are given in Tables (1) and (2) respectively. 

Figure (5-a) shows the time history of the amplitudes of transverse 
vibrations of the steel cantilever beam when controlled by the MIKSC 
method using one PZT piezo-electric actuator, bonded to the element near 
the fixed end of the beam. In this figure the effect of bonding layer on 
the physioal properties of the system is neglected. The performance of 
the system under such conditions will serve as a datum for measuring the 
actual effect of the bonding layer when it is accounted for in the 
system's model. 

With the time sharing concept, the HIMSC utilizes effectively the 
installed aotuator such that it provides the necessary action to oontrol 
the dominant modes and then shares the controller among the other modes 
until the vibrations of the system is completely damped out based on the 
maximum modal energy ranking. This time sharing concept can best be 
understood by considering Figures (5-b) and (5-c). 

Figure (5-b) shows the control mode that has the highest modal 
energy at any instant of time. It can be seen that the actuator is used 
first to attenuate the amplitude of vibration of the first mode, which 
currently has the highest modal energy. After a small time interval the 
control action is switched to control the second mode since its modal 
energy becomes higher than the energy of the other five modes. This 
action of time sharing the single actuator between the modes continues 
until all modes are brought under control. This is demonstrated clearly 
in Figure (5-c) by the continuously decaying vibration energy of the 


system. 



Table (1) - Properties of bonding layers [10] 


Bonding Layer 

! Density (Kg/m^) | 

Young’s Modulus (MN/m 2 ) 

laotac 1 

1 890 I 

1.1 

Eastman 910 2 

1 1050 | 

1780.0 


1 Rubber adhesive tape nade by Minnesota Mining and 
Manufacturing. St. Paul, Minn. 

2 Cyanoacrylate adhesive made by eastman Kodak Co., Rochester, 
N.T. 





Table (2) - Properties of a typical piezo-electric actuator [8] 


Material 

1 

PZT 

-Che.. <?e coefficient d 
10“ (m/v) 

1 

1 

123 

Young's modulus 
(G N/m 2 ) 

1 

1 

139 

Max. voltage v* 
(Mv/m) 

1 

1 

1 

Max. tensile strength o* 
(M N/m 2 ) 

1 

1 

45 

Density 

(kg/m*) 

1 

1 

7500 








Figure (5-a) - Tine response of steel cantilever bean controlled by one 
PZT actuator negleoting effect of bonding 



Figure (5-b) - Control node of highest nodal energy for the cantilever bean 



Figure (5-c) - Total nodal energy of oantllever bean 





(A) EFFECT OF BORDUG LATER TBICOESS 


Figures (6-a), (6— to) and (6-c) display the time histories of the 
amplitudes of transverse vibrations of the cantilever beam when one PZT 
piezo-electric actuator is bonded to element 3 using an Isotac bonding 
layer of thickness of 0.42mm, 0.84mm and 2.1mm respectively. In 
dimensionless form, the selected bonding layers have thickness equal to 
0.2, 0.4 and 1.0 that of the actuator thickness respectively. 

A comparison between Figures (5-a) and (6-a) indicates that the 
maximum amplitude of vibration has decreased from 1.43E-5m to 1.39E-5m 
when the elastic and Inertial properties of the bonding layer, that has a 
thickness of 0.00042m, are accounted for. Further decrease in the maximum 
amplitude of vibration is observed as the thiokness of the bonding layer 
is increased as displayed in Figures (6-b) and (6-o). 

Detailed analysis of the figures is summarized in Table (3) based on 
the displacement and control force indices U d and 0 C which are defined by 
the following expressions: 


°d 8 



t=0 i=1 
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2 2 F4 2# At 

t=0 i=1 


( 16 ) 


(17) 


where 


N is the number of d.o.f. of system 
at is the integration time increment 
t* is the maximum time limit of integration 


Table (3) indicates clearly that including the effect of bonding 
layer results in decreasing the displacement index U d monotonely but on 








Table (3) - Effect of bonding layer thickness on the displacement 
and control force Indices for a steel cantilever beam 
controlled by one PZT actuator at element 3 


Thickness I Displace^r*: index I Control force index 
(m) I (x10 lu ) I (xIO 7 ) 


0 

1 

2.8244 

1 

1.8128 

0.00042 

1 

2.6505 

1 

3.0164 

0.00084 

1 

2.5531 

1 

3.4443 

0.00210 

1 

2.3787 

1 

9.6501 


the expense of a significant increase in the required oontrol voltages. 
For example, using bonding layer thickness of 0.00042m reduces the 
displacement index from 2.8244E-10 to 2.6505E-10 whereas it increases the 
control force index from 1.8128E-7 to 3*0l64E-7, which amount to an 
increase of 66%, 

The reduction in the amplitudes of vibration of the beam resulting 
from using thicker bonding layers is attributed primarily to the 
stiffening effect produoed by the bonding layers to the original 
actuator-beam system. However, using such thicker bonding layers would 
require large control forces to damp out the vibrations of the flexible 
beam. This is reflected directly into high control voltages as can be 
seen clearly from Figures (7-a)» (7-b) and (7-c) for Isotao bonding layer 
thickness of 0.42mm, 0.84mm and 2.1mm respectively. 

Therefore, we have a push and a pull situation where thicker bonding 
layer are preferred to minimize the amplitudes of oscillation and thinner 
bonding layers are desirable if one is to limit the magnitude of the 
control voltages. 

(B) EFFECT OF BOND DIG LITER MATERIAL 

The effect of changing the material of the bonding layer from Isotac 
to Eastman 910 on the displacement and control force indices is 
summarized in Table (4). The results listed in the table are obtained 
when one PZT actuator is bonded to element 3 and a 0.42mm thick bonding 
layer is used. 

Table (4) indicates that changing the bonding layer material from 
the soft Isotac to the more stiff and dense layer of Eastman 910 results 
in insignificant reduction in the displacement index but in a 
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Time history of voltage controlling a steel cantilever 
beam by one PZT actuator bonded to element 3 with Isotac 


of different thickness 




Table (4) - Effect of bonding layer material on the displacement 
and control force indices for a steel cantilever beam 
controlled by one PZT actuator at element 3 


Material 


| Displacement index - ! 
I (xIO 10 ) I 


Control force index 
(xIO 7 ) 


2.6505 

2.6499 


considerable increase , of about 20. 4f, in the forces necessary to 
oontrol the vibration of the beam. 

(C) EFFECT OF ACTUATOR LOCATION 

Figures (8-a) and (8-b) show the time historlES of the amplitudes of 

transverse vibrations of the steel cantilever beam when a PZT actuator is 

bonded to element 2 and 3 respectively. The bonding layer considered, in 
this case, is Isotac vhioh has a thickness of 0.00042m. 

The associated displacement and control force indices are listed in 
Table (5). 

Table (5) as well as Figures (5-a), (8-a) and (8-b) indicate that it 

is preferable to plaoe the actuator at element 3 near the fixed end of 

the beam in order to obtain the minimum displacement index. On the other 
bead, if the minimum oontrol forces are concerned, then this actuator 
should bonded to element 1. 


COBCLUSIOIS 

This paper has presented a comprehensive analysis of the control of 
vibration of simple structural elements using piezo-electric actuators. A 
Modified Independent Modal Space Control (MIMSC) method is utilized to 
actively control the vibration of the flexible system. 

The effect of the physical parameters of the actuator and its 
bonding layer on the elastic and inertial characteristics of the flexible 
beam are included in the analysis. It is observed that accounting for the 
elastic and inertial properties of the bonding layer results in reducing 





Table (5) - Effeot of actuator looation on the displacement and 
oontrol force lndloes for a steel cantilever beam 
controlled with one PZT actuator 


Actuator 

location 

1 

1 

Displacement index 

1 

1 

Control force 
(x10 7 ) 

Element 1 

1 

5.2855 

1 

0.1021 

Element 2 

I 

3.3681 

1 

2.5494 

Element 3 

1 

2.6505 

1 

3.0164 


the amplitudes of vibrations of the structural element but on the expense 
of higher oontrol forces. 

The study has demonstrated also the effectiveness of the HIMSC 
method in controlling the vibration of flexible beams that have large 
number of degrees of freedom with very small number of actuators. 
Furthermore, the presented method can be readily extended to large 
structures without any modifications. 
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